Noname manuscript No. 

(will be inserted by the editor) 



Coupled-resonator optical waveguides: Q-factor and 
disorder influence 

Jure Grgic, Enrico Campaioli, S0ren Raza, 
Paolo Bassi and Niels Asger Mortensen 



f****) Received: date / Accepted: date 

Abstract Coupled resonator optical waveguides (CROW) can significantly reduce 
I | light propagation pulse velocity due to pronounced dispersion properties. A number 

. of interesting applications have been proposed to benefit from such slow-light propa- 

^f*j ' gation. Unfortunately, the inevitable presence of disorder, imperfections, and a finite 

Q value may heavily affect the otherwise attractive properties of CROWs. We show 
^ ' how finite a Q factor limits the maximum attainable group delay time; the group index 

O . is limited by Q, but equally important the feasible device length is itself also limited 

by damping resulting from a finite Q. Adding the additional effects of disorder to this 
picture, limitations become even more severe due to destructive interference phenom- 
ena, eventually in the form of Anderson localization. Simple analytical considerations 
' demonstrate that the maximum attainable delay time in CROWs is limited by the 
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intrinsic photon lifetime of a single resonator. 



1 Introduction 



Slowing down the speed of light enhance phenomena like nonlinearities [T], gain/ab- 
sorption [2], [3], and phase sensitivity [4]. It can be also useful for practical applications 
like delay lines, optical memories, and low threshold lasers [5J. CROW structures offer 
particular ways of slow light guiding, with photons hopping sequentially from one 
, resonator to the next. The concept and proposal of CROWs were introduced in 1999 

by Yariv et al. [6]. At that point it was obvious that such structures would yield big 
interest in the optical community due to the potential application possibilities 7,8,9, 
110] . In particular, it is clear that near the band edge the group velocity is significantly 
reduced in the ideal and lossless structure. Unfortunately, a certain level of damping 
will always be present due to radiation losses, intrinsic losses of materials, and any 
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other mechanism that may dissipate or scatter part of the electromagnetic energy. Any 
fabrication processes will introduce some variations in the properties of the individual 
resonators that will serve to dissipate and scatter light propagating in the CROW. 
Such non-uniformities can indeed seriously affect slow light propagation jlllll2l l71ll3|. 
emphasizing the importance of quantifying their effect and unavoidable consequences. 

When discussing the potential of slow-light waveguides the emphasis is often on the 
attainable group index n g — c/v g or the group velocity v g . However, in applications 
involving delay lines and buffers, the group delay r = L/v g is the key parameter. 
Obviously, the longer the waveguide the longer the group delay! However, this trivial 
statement implicitly neglects the decay of the pulse as it propagates down the waveguide 
and eventually the pulse has lost its initial strength and intensity. In this paper, we 
treat the issues of delay and decay on an equal footing by emphasizing that the decay 
length £ serves as an effective cut off for L (section 2), leaving us with a maximal 
attainable group delay r max of the order £/v g [X41IT5] . 

The manuscript is organized as follows. From coupled- mode theory (section 3), 
our key observation is that the maximum attainable delay time in CROWs is limited 
by the intrinsic photon lifetime of a single resonator (section 3.1). The presence of 
disorder and scattering will further serve to reduce this bound (section 3.2), as shown 
numerically with the aid of a Green function method (details given in appendix). 
Finally, conclusions are given (section 4). 



2 Delay versus decay 

To facilitate quantitative predictions of Tmax we imagine a CROW of length L and the 
transmission T through this segment is then conveniently parameterized by 

TH = exp(-iM). (1) 

Obviously, this parametrization excellently represents the exponential decay of the 
power associated with absorption and other loss mechanisms captured by the finite Q 
factors of the resonators. Likewise, in the presence of localization, the average trans- 
mission also has an exponential distribution with the scale given by the localization 
length. Thus, the length scale — —L\nT{ui) captures the combined effects of 

disorder-induced localization and other loss accounted for by the finite Q-factor. 

In the context of pulse delay, £ represents an estimate of the maximal length of 
the CROW that we could imagine in any practical application. Extending the length 
L of CROW beyond £ would effectively suppress the output power, thus jeopardizing 
any benefits of slowing down a wave package. The maximal group delay is thus a 
balance between a slow group velocity and a long propagation length, i.e. r max (uj) = 
£,(uj)/v g (oj). Since the group velocity is also inversely proportional to the density-of- 
states we may conveniently rewrite this expression as 

Tmax(w) = — Tvp(uj) lnT(oj). (2) 

The key result of this paper is that Tmax is always limited by the single-resonator 
photon lifetime, i.e. 

T max (oj) < t p (3) 
where the equality applies to the situation with absence of disorder. 
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In the following we outline our account for this result by coupled-mode theory 
and in particular, we study first the limit of ideal CROWs (but with finite Q) and 
subsequently we address the effects of disorder. 



3 Coupled-mode theory 

We consider a chain of coupled optical resonators, where the jth resonator, if left 
isolated from the other resonators, is characterized by a resonant field 

Ej (r, t) = Ej (r) exp [i{Q 3 + iTjfflt] , (4) 
where f2j is the resonance frequency and Fj represents the resonance line width. The 
energy in the resonator Ej (uj) then has a Lorentzian frequency distribution corre- 
sponding to the density-of-states 

PjH = i(a J -^ + (r,/2)2- ^ 

The associated quality factor Qj = fij/Tj may conveniently be parametrized as a 
photon life time t p = Qj fi. 

Next, imagine a chain of coupled resonators of the above kind. We follow the work 
of Yariv and co-workers [6] and write the electrical field as a linear combination of 
the isolated resonator fields, i.e. Ej(r) = ^ • ipjEj(r), where the expansion coefficients 
are denoted by tp to emphasize the similarities with the notation in Ref. [16III 7] for 
the associated problem of electrons in a quantum wire. In the regime of weak nearest- 
neighbor coupling the equations linearize to 

(fij + Wj/2)ipj - 7j+i jtp J + 1 - lj-l.jipj-l = wV>j (6) 

which has a form resembling the tight-binding chain in condensed matter physics |18l 
I16| . To further emphasize this connection we write the coupled equations in a matrix 
form, i.e. Hip = ujip, where the 'Hamiltonian' matrix H has elements 

fly = {Qj + iTjfflSy - JljSi ±1)j (7) 

with 7y = 7*; so that the off-diagonal part of H is Hermitian. 



3.1 The influence of finite Q factor 

In the case where the resonators are all identical and arranged in a fully periodic se- 
quence, we may without loss of generality suppress all indices, thus making further 
analytical progress possible. The electromagnetic states now form a continuous fre- 
quency band and the problem is easily diagonalized by the Ansatz ipj+1 — exp(ina)tpj . 
The resulting dispersion relation is of the form 

w(k) = Q (\ + i ^Q^j - 27cos(fta) (8) 

where k = k' + ik" is the complex valued Bloch wave vector and a is the lattice constant 
of the periodic arrangement of resonators. Eq. ((SJ) corresponds to the theory by Yariv et 
al. [B] to also include resonators with a finite Q-factor|j 

1 To ease the comparison to our previous work in Ref. 1141 we note that lu(k) = 
Q (l + *2^) [1 — 7cos(rea)] where 7 is for simplicity considered real and given by 7 = 2-y/Q. 
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The group velocity may formally be calculated from the dispersion relation in 
Eq. |8]l. We imagine the situation where the CROW is excited by a monochromatic 
laser with a well-defined frequency of the light. Care should thus be taken that uj is 
to be considered real while k, may be complex, and the group velocity is then formally 
given by 



Along the same lines, we may also calculate the density-of-states from the dispersion 
relation in Eq. (JH]). We emphasize, that for the particular case of a one-dimensional 
chain, the density-of-states is inversely proportional to the group velocity, 



Clearly, a vanishing group velocity will be associated with a diverging density-of- 
states and vice versa. This also explains why a light pulse can not completely come to 
stop in a real system, since no real system will exhibit a true singular density-of-states 
the inevitable presence of absorption, radiation, and imperfections will serve as effective 
broadening mechanisms. Nevertheless, to illustrate the basic physics, it is common to 
consider lossless structures and only focus on the real part of the dispersion properties, 
since the damping is anyway assumed to be modest. However, in the present context of 
slow-light propagation the two issues are not easily separable and care must be taken. 

Figure [T] illustrates the contrast between a CROW made from lossless resonators 
(dashed lines) and finite-Q resonators (solid lines), respectively. The left panel illus- 
trates the usual cosine-band dispersion relation, i.e. the relation between the frequency 
10 and the real part n' of the complex- valued Bloch wave vector k, = k + i/t . Likewise, 
the middle panel illustrates the corresponding damping, i.e. the relation between fre- 
quency u) and the imaginary part k" of the Bloch wave vector. Finally, the right panel 
shows the density-of-states associated with the dispersion diagram in the left panel. 
The difference between the ideal structure (Q — > oo) and one employing resonators of 
finite Q is easily contrasted by comparing the dashed and solid lines, respectively. We 
emphasize that the main cause of a finite Q factor is to smear out van Hove singular- 
ities in the density-of-states and to weaken the slow down of light pulses propagating 
near the band-edges of the CROW. Of course, the finite Q also introduces damping 
throughout the entire band, though be most pronounced near the band edges due to 
slow-light enhanced absorption [19]. Finally, we note that quite steep bands appear 
outside the traditional band of extended states, though of course with a significant 
attenuation as evident from the middle plot illustrating the k" dependence. 

In our previous work we have carefully discussed the influence of the finite Q factor 
on the saturation of the group index [1411 15] . The main conclusion is that in the center 
of the band the group velocity is (to lowest order in 1/Q) insensitive to the finite photon 
life time. On the other hand, at the band edges, initially supporting pronounced slow 
down, the group velocity scales quite unfavorably with the Q factor, making the slow- 
light regime challenging to explore. 

Here, we focus on the attainable group delay when the slow down and the damping 
is treated on an equal footing. For the ideal CROW, the group delay is given by 
r = L/vg . However, in the presence of a finite Q, the length L is effectively cut off by 
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the damping length £ = 1/2k" associated with the exponential decay in Eq. ([TJ) . This 
gives an upper bound and Eq. ([2} may in this case formally be rewritten as 

r max H = ^7^. (11) 
Combining the full results for k! and k" and expanding in 1/Q we get [14] 

TmaxH =T P + 0(Q- 1 ). (12) 

The main conclusion from this analytical exercise is that the maximal group delay is 
limited by the photon life time t p = Q/£2 of the individual resonators, independently 
on the actual frequency. While being a quite intuitive results, it has important and 
overlooked consequences with respect to how much delay one may envision in future 
designs of optical buffers and delay-line architectures. Despite the reduced group veloc- 
ity near the band edges, the advantage of a slowly advancing wave package is balanced 
by a reduced propagation length, see the middle panel of Fig. [1] We emphasize that 
compared to a single resonator, the CROW may of course offer the advantage of an 
increased bandwidth. Likewise, the strongly suppressed group-velocity dispersion at 
the band center might also be beneficial in some applications. 

In the following we discuss how this result is modified in the presence of disorder. 
However, we may at this stage anticipate that disorder may only further limit time 
delay, thus in general adding the 'lesser sign' in Eq. © as compared to the equality 
derived in Eq. (JT2J ■ 



3.2 The influence of disorder 

We next turn to disordered waveguides, formally allowing for CROWs composed of res- 
onators with a resonator-to-resonator fluctuation in the resonance frequency Qj and 
the linewidth Fj as well as fluctuations in couplings 7y between neighboring resonators. 
For simplicity we will neglect fluctuations in the linewidth so that all resonators have 
the same Q factor. For the resonance frequencies and the couplings we will further 
assume uncorrelated Gaussian distributions P(Qj) and P{"yij), respectively. The dis- 
tributions have mean values f2 and 7 corresponding to the ideal CROW parameters, 
while fluctuations cjq = VS 2 f} and 07 = \J 5 2r y around the mean values can be varied 
to mimic different strengths of disorder. 

We employ a Green's function method to calculate the transmission and the density- 
of-states which allows us to evaluate Eq. (J2J) for each member of the ensemble. The 
Green's function is formally obtained by inverting the matrix H. However, since its 
dimensions are formally infinite we imagine a segment of disordered CROW (containing 
N resonators) sandwiched between to semi-infinite ideal CROWs, i.e. with no disorder. 
With the aid of Dyson's equation this apparently unsolvable problem can fortunately 
be turned into a matrix problem of finite dimension (corresponding to the dimension 
N of the disordered segment) and the retarded Green functions of the CROW can now 
be found from a numerical inversion of a sparse N x N matrix problem [171116] 

g(u)=[uI-H-E(u)]~ 1 , (13) 

where J is a unit matrix and the couplings to the two semi-infinite ideal CROWs are 
accounted for by a complex-valued frequency dependent self-energy S. The details of 
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this approach are given in the appendix [X] which lists expressions how to obtain the 
transmission T and the density-of-states p from the Green's function Q. 

Once the Green function is obtained, the maximal time delay, Eq. @, is thus easily 
evaluated with the aid of Eqs. (|16|l and (|18|) . We take advantage of standard matrix 
inversion routines to numerically study the statistical properties of large ensembles of 
disordered CROWs. In principle this allows us to study statistical moments to any 
order, but for simplicity we will here focus on the average properties (first moment) 
and only highlight the CROW-to-CROW fluctuations (second moment) by displaying 
results for particular members, chosen randomly from the large CROW ensemble. 

In the panel b) of the Figure [2] we show results for the ensemble-averaged DOS 
(blue lines). It is clearly seen how disorder, in addition to a finite Q, serves to further 
broaden the ensemble-averaged DOS near the band edges. Comparing these results to 
the DOS associated with one particular realization of the disordered CROW (red line), 
it is however clear that pronounced sample-to-sample fluctuations are to be expected. 
In particular, the formation of Anderson localized states near the band edges, due to 
the strong interference of light waves, is apparent. These fluctuations in the DOS are 
quite naturally inherited by other central quantities, such as the maximal group delay 
Tmax and the transmission T. In panel c) the value of r m ax is normalized with photon 
lifetime in single resonator t p . For the ensemble-averaged maximal delay time, disorder 
is seen to further suppress r ma x below the bounds by r p . However, from a practical point 
of view it is alarming to see fluctuation comparable to the mean value, as indicated by 
the strongly fluctuating results for a particular realization of the disorder (red line). 
For the transmission in panel d) we see a similar picture with a strong suppression of 
the transmission near band edges, but with the pronounced transmission fluctuations 
appearing throughout the entire band. For comparison, the dashed line shows the result 
of unity transmission (T = 1) for an ideal CROW with infinite Q, while the green line 
shows the pronounced suppression of transmission in the presence of a finite Q, but in 
the absence of any additional disorder. The quite abrupt drop in transmission near the 
band edges is associated with slow-light enhanced absorption [20] as compared to the 
center of the band where slow-light enhancement is almost absent. 

There is an interesting interplay between the finite Q factor and the amount of 
disorder in the structure. In Figure[3]we plot T m ax/T p as a function of disorder strength 
a I Q = <t 7 /Q = oq /!?, evaluated for a frequency corresponding to the band center. The 
different curves represent different values of the Q factor. Quite intuitively, the lower 
the Q value the less sensitive is the result to disorder, keeping in mind that Tmax, T, 
and n g themselves would be heavily suppressed in the presence of a low Q. For higher 
Q values, it consequently becomes increasingly challenging, in terms of disorder, to 
take full advantage of the high intrinsic photon life time Tp. 

4 Conclusion 

In conclusion, we have derived an explicit relation for the dispersion relation of CROWs 
made from resonators with a finite Q factor. A finite Q profoundly influences the 
van Hove singularities near the band edges with a resulting limitation of the group 
index while at the center of the band the dispersion properties are less affected. Simple 
analytical expressions are supported by calculations of the group velocity, demonstrat- 
ing how the Q enters on an equal footing with the coupling 7 corresponding to the 
competing time scales associated with photon decay and tunneling. In the context of 
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practical applications involving the group delay, we note that the maximal attainable 
group delay appears as a balance between the reduced group velocity and the the decay 
length. Explicit calculations show that irrespectively of the underlying bandstructure, 
the maximal group delay is limited by the photon life time of the resonators. This 
illustrates the importance of addressing propagation loss and slow-light on an equal 
footing. Any inevitable presence of disorder will serve to further suppress the attainable 
group delay and pronounced sample-to-sample fluctuations may arise. 
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A Details of Green's function approach 

The self energy in Eq. (1 1 3 P is given by S = S L + S R where the contributions from the left 
and right semi-infinite CROWS are given by 

{ S p( E )}jl = -7cxp [t«(w)o] S j3p S Spl (14) 

with sl = 1 and = N. The wave vector is related to the energy through the usual cosine 
dispersion relation derived above, see Eq. JSJ, corresponding to 



exp [i K {uj)a] = — 4 a ■ (15) 

The transmission probability may now conveniently be written as a trace formula 

7V) = Tr[rj,(u>)g(u)r R (u)gi(u)], (w) 

where 

r r (u) = i[S p (u)-l!l(u)]. (17) 
Likewise, the total density-of-states (per resonator) is given by 

1 N 

P H = -^ft-H (18) 
i=i 

with the corresponding local density-of-states governed by the diagonal part of the Green's 
function, 

PM =--M£?jj(w)}. (19) 
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Fig. 1 Complex dispersion relation for a CROW. Dashed lines are for Q = oo while solid 
lines correspond to Q = 10 2 . The left panel shows the frequency w versus the real part of the 
Bloch wave vector re', the middle panel shows the frequency uj versus the imaginary part of 
the Bloch wave vector re", and the right panel shows the density-of-states p (per resonator). 
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Fig. 2 Properties of a disordered CROW, with blue lines indicating ensemble- averaged prop- 
erties while the red lines illustrate the results for a particular realization of the disorder, thus 
emphasizing pronounced CROW-to-CROW fluctuations. Panel (a) shows the local DOS pj for 
a particular realization of the disorder and panel (b) shows the corresponding results for the 
total DOS p (per resonator). Panel (c) shows the maximal group delay r max . Panel (d) shows 
results for the transmission. For comparison, the dashed line shows the unity transmission for 
an ideal crow, while the green line is for a non-disordered CROW, but with a finite Q. 



12 




